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SECTION  1 


INTRODUCTION 

During  the  past  decade  the  field-effect  transistor  has 
become  the  single  most  important  microwave  solid-state  device.  It  is 
used  in  small-  and  large-signal  amplifiers,  low-noise  amplifiers, 
oscillators,  mixers,  and  a  variety  of  control  applications.  The  FET 
has  been  the  subject  of  numerous  papers,  leading  to  a  better  under¬ 
standing  of  its  operation  and  potential.  Early  FET  papers  were 
analytic  approximations  which  neglected  the  effect  of  space  charge 
in  the  channel.  With  faster  computers  and  better  numerical  methods, 
full  solutions  of  the  two-dimensional  carrier  transport  and  Poisson's 
equations  became  possible.  Results  from  these  finite-difference  and 
finite-element  solutions  increased  our  understanding  of  the  device. 
However,  the  microwave  FET  structure  poses  a  basic  problem  for  the 
device  modeler.  An  X-band  FET  has  a  typical  doping  density  of  1017/cm3, 
with  higher  doping  for  higher  frequency  operation.  The  characteris¬ 
tic  length  in  a  semiconductor  material  is  the  debye  length  which  is 
given  by1 


where  k  =  the  Boltzmann  constant, 

T  =  the  temperature, 
e  =  the  dielectric  constant, 
q  =  the  electronic  charge, 

N,  =  the  material  doping 
a 


(i) 


i 


and  the  characteristic  time  is  the  dielectric  relaxation  time  which 


is  given  by 


£  1 

q  pNd  * 


(2) 


where  \i  =  the  low  field  mobility. 

For  of  1017/cm3,  the  debye  length  is  approximately  0.02  pm 
and  the  dielectric  relaxation  time  is  approximately  0.01  ps.  The  prob¬ 
lem  with  dc  modeling  in  two  dimensions  is  constrained  by  the 
small  mesh  size  and  the  resulting  huge  number  of  calculations  required 
to  get  a  solution.  The  RF  or  switching  problem  means  doing  a  solution 
approximately  100  times  per  ps.  Of  course,  though  various  approxima¬ 
tions  and  clever  programming,  these  problems  can  be  overcome  to  some 
extent. 

The  approach  used  in  this  report  is  a  combination  of  the  ana¬ 
lytic  and  the  full  simulation.  The  Poisson  and  electron  transport  equa¬ 
tions,  altered  to  account  for  the  two-dimensional  nature  of  the  problem, 
are  solved  in  one  dimension.  Several  analytic  expressions  are  used  to 
determine  the  geometry  of  the  channel.  These  approximations  give  a  set 
of  equations  that  can  be  solved  rapidly.  At  the  same  time  most  of 
the  two-dimensional  device  physics  remains. 

The  approximations  used  in  the  FET  model  and  the  resulting 
equations  are  given  in  Section  2.  The  model  is  used  to  describe  FET 
operation  and  to  compare  it  with  a  full  two-dimensional  simulation  in 
Section  3.  Additional  studies  of  GaAs  FETs  with  varying  velocity-field 
and  diffusion-field  characteristics  are  given  in  Section  k,  where  the 
effect  of  gate  length  is  also  studied.  A  summary  of  the  major  results 
is  given  in  Section  5- 
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SECTION  2 


MODEL  DESCRIPTION 

The  purpose  of  this  work  is  to  devise  a  simplified  but  realistic 
model  for  the  study  of  field-effect  transistors.  There  are  a  variety 
of  approaches  to  FET  modeling.  Shockley’s2  original  paper  used  the 
gradual  channel  approximation  and  the  approximation  of  zero  channel 
space  charge.  This  type  of  analytic  model  was  also  used  by  Pucel 
et  al.,3  Shur  and  Eastman,1*  and  Dawson  and  Frey5  to  further  study  the 
FET.  These  models  all  neglect  space  charge  and  other  physical  effects 
within  the  conducting  channel.  At  the  other  extreme  of  complexity 
is  the  full  two-dimensional  solution  of  the  carrier  transport  and 
Poisson's  equations.6-10  The  model  developed  here  is  a  compromise 
between  the  two  types.  The  hope  is  to  have  a  model  that  is  simple 
(i.e.,  low  cost)  enough  to  run  a  variety  of  cases  and  still  complex 
enough  to  provide  a  further  understanding  of  the  properties  and  limita¬ 
tions  of  FETs.  The  model  is  intended  to  provide  most  of  the  information 
which  requires  full  two-dimensional  simulation. 

The  approach  used  can  be  discussed  with  the  help  of  Figure  1. 

This  shows  the  voltage  distribution  and  gate  depletion  layer  geometry 
for  three  FET  structures  obtained  with  a  two-dimensional  simulation. 
Throughout  the  conducting  channel,  the  voltage  contours  are  approximately 
straight  and  perpendicular  to  the  axis  of  the  device.  This  is  true  in 
the  limits  of  long  gates  and  is  the  basis  of  Shockley's2  gradual 
channel  approximation.  However,  as  can  be  seen  from  Figure  1,  the 
approximation  of  equipotential  lines  perpendicular  to  the  axis  of  the 


-3- 


I 


i 


FIG.  1  CALCULATED  VOLTAGE  DISTRIBUTION  WHEN  V  =  0  an 

sg 

V  .  =  5.0  V  (FIELD-DEPENDENT  ELECTRON  MOBILITY) 
sa 

(a)  2W/W  =  5.0,  (b)  2W/W  =  1.0,  (c)  2W  /W 
6  c  &  c  &  ^ 

(KENNEDY  AND  O'BRIEN*) 


device  is  also  reasonable  for  very  short  gate  devices.  This  approximation 
is  equivalent  to  assuming  that  there  is  no  vertical  electric  field  in 
the  channel.  With  the  additional  assumption  of  an  abrupt  transition 


between  the  conducting  channel  and  the  depletion  layer  a  set  of  differ¬ 
ential  equations  can  be  written  to  describe  the  channel. 

Poisson's  equation  and  the  electron  continuity  equation  for  the 
channel  can  be  written  using  the  geometry  and  notation  of  Figure  2a. 

The  electric  field,  voltage,  electron  concentration,  and  channel  geom¬ 
etry  factor  B  are  functions  of  position  x  along  the  channel.  The 
velocity  and  diffusion  are  functions  of  position  through  the  velocity- 
electric  field  and  diffusion-electric  field  curves  for  the  device 
material.  B(x),  the  channel  geometry  factor,  is  normalized  to  the 
height  of  the  channel  a.  A  small  volume  element  of  length  Ax  and  unit 
height  is  shown  in  Figure  2b.  If  it  is  assumed  that  there  is  (a)  no 
current  flow  at  the  bottom  (this  is  the  channel-substrate  interface), 
(b)  no  current  flow  at  the  sloping  top  (this  is  the  channel-depletion 
layer  interface),  (c)  an  average  electron  concentration  within  the 
element  of  and  (d)  an  average  electric  field  on  the  right  and  left 
sides  of  Ej+1  and  E^,  then  Poisson's  two-dimensional  equation 

V.E  =  -  (N  -  Nj  (3) 

e  a 

can  be  written  as 


biei  - 


bi+iei+i 


f  (Ni  -  — 'ir ~ 


ik) 


For  small  Ax  this  is  approximately 
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GATE  SCHOTTKY 


BI  BI+I 

(b) 

FIG.  2  (a)  FET  GEOMETRY  AND  NOTATION,  (b)  VOLUME  ELEMENT 


GEOMETRY  AND  NOTATION. 


3 (BE)/3x 


(5) 


^  (N  -  N  )B  . 
e  a 


Similarly ,  the  electron  continuity  equation  becomes 


Vi  -  V1V1 


ANI  AXBI+1  +  BI 


where 


3 (BJ)  _  3 (BN) 

3x  "  3t 

J  =  the  electron  current  density 

=  Qv[E(x)]N(x)  -  QD[E(x)][dN(x)/dx]  , 


v[E(x)]  =  the  velocity  and 
D[ E(x ) ]  =  the  diffusion  coefficient. 
The  voltage  along  the  channel  is 


•x 

J  x=c 


E(x'  )  dx'  . 


The  channel  opening  B(x)  depends  on  the  voltage  along  the  channel. 

The  current  in  the  channel  depends  on  the  material  parameters 
used.  The  form  of  the  velocity-field  and  diffusion-field  curves  are 
from  Bauhahn.11  The  curves  are 


1  +  (E/E  )** 


1  ♦  (E/Ed  >“ 


where  =  the  low  field  mobility, 

v  =  the  saturated  velocity, 
sat 


! 


E^  =  the  peak  field  for  velocity, 

Dsat  =  the  his*1  field  diffusion  coefficient  and 

E,  =  the  corner  field  for  diffusion, 
d 

Before  Equations  5  and  7  can  be  solved,  an  equation  relating  the 
voltage  along  the  channel  and  the  depletion-layer  width  or  channel 
geometry  factor  B  is  needed.  The  voltage  drop  across  the  depletion 
layer  is  related  to  its  width  by 

Vd  =  (q/2e)Ndu)2  .  (12) 

One  possible  approximation  is  to  assume  that  the  equipotential  lines 
within  the  depletion  layer  are  parallel  to  the  gate  of  the  Schottky 
barrier,  as  shown  in  Figure  3.  This  approximation  is  typically  used  in 
analytic  FET  models  and  is  correct  in  the  limit  of  long  gates.  With 
this  approximation,  the  voltage  across  the  depletion  layer  depends  on 
the  voltage  difference  between  the  gate  and  the  channel.  Equation  12 
can  be  used  to  find  the  depletion-layer  width  and  the  resulting  channel 
geometry  factor  B(x).  This  approximation  is  referred  to  as  the 
parallel-potential  approximation. 

The  problem  with  the  parallel-potential  approximation  is  that  it 
always  overestimates  the  width  of  the  gate  depletion  layer  and  under¬ 
estimates  the  channel  geometry  factor  B(x).  This  can  be  shown  by 
referring  to  Figure  1.  Although  the  constant  voltage  lines  on  the 
conducting  channel  are  approximately  parallel,  they  bunch  together  in 
the  depletion  layer  near  the  drain  end  of  the  gate  depletion  layer. 

Thus  a  sloping  straight-line  approximation  can  be  made  for  the  equi¬ 
potential  lines  within  the  depletion  layer.  If  the  voltage  along  the 
top  of  the  channel  outside  the  depletion  layer  is  known 
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(from  Equation  9),  Equation  12  can  be  used  to  find  the  corresponding 
depletion-layer  width  between  the  gate  and  the  drain.  Although  the 
geometry  of  the  equipotential  lines  is  a  two-dimensional  problem  within 
the  depletion  layer,  there  is  only  an  x-direction  electric  field  along 
the  top  surface  of  the  channel.  The  voltage  along  the  surface  is 


V*  (x)  = 


V  +  AV.  +  (Q/2e)N.Ax2 


r»  _  Sl-  V  •  2 

2e  d 


I 


f 


and  the  depletion-layer  width  at  x  is 

ui'  (x)  =  r  sin  6  .  (l6) 

Using  Equations  15  and  l6  yields  the  depletion-layer  width: 

u)'(x)  =  o)(x)  sin  0/0  .  (17) 

This  set  of  equations  can  be  used  to  find  the  electron  and  field 
distribution  along  the  channel.  The  program  starts  from  an  initial  guess, 
typically  with  n(x)  equal  to  the  background  doping  and  B(x)  =  1.  For 
a  given  gate  and  drain  voltage,  Equation  5  can  be  solved  to  find  E(x). 

Then  Equation  9  can  be  integrated  to  find  the  voltage  along  the  channel. 
The  channel  geometry  can  be  found  using  the  voltage.  When  the  geometry 
and  electric  field  are  known.  Equation  7  can  be  used  to  update  the 
electron  distribution.  This  process  is  continued  until  a  convergent 
solution  is  obtained. 

Two  versions  of  this  program  have  been  programmed;  one  in  the 
HPL  programming  language  running  on  an  HP  9825  calculator  and  one 
in  FORTRAN  running  on  an  AMDAHL  computer.  The  HPL  version  requires 
several  minutes  to  find  a  converged  solution.  A  complete  set  of  solu¬ 
tions,  over  a  range  of  drain-source  and  gate-source  voltages,  requires 
less  than  a  second  of  AMDAHL  cpu  time. 


SECTION  3 


COMPARISON  OF  THE  SIMPLE  MODEL  AND  THE  TWO-DIMENSIONAL  RESULTS 

Introduction .  The  operation  of  field-effect  transistors  is 
described  in  this  section.  A  general  discussion  of  FET  operation  is 
given  first  so  that  the  factors  effecting  the  device  performance 
can  be  better  understood.  Then  the  results  of  the  present  work  are 
compared  with  owo-dimensional  results.  This  will  show  the  approxi¬ 
mations  and  accuracy  of  the  present  model. 

FET  Operation.  The  operation  of  FETs  can  be  explained  using 
Equations  5  and  7.  If  the  electron  concentration  in  the  channel  is 
equal  to  the  channel  doping  and,  in  the  case  of  dc  operation,  the 
equations  become 

B(x)E(x)  =  constant  (l8) 

and 

B(x)j(x)  =  J  ,  (19) 

dc 

the  constant  in  Equation  l8  is  set  by  the  applied  drain-source  voltage. 
The  current  density  is 

J(x)  =  qNduQ(x)E(x)  .  (20) 

The  current  density  field,  and  geometry  are  shown  in  Figure  5  for 
the  simple  approximation  of  a  rectangular  depletion  layer.  The 
voltage  across  the  device  is 

V  =  Ejtxj  +  (L  -  x2)  +  (x2  -  Xj/b2)]  .  (21) 
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The  current  is  related  to  the  field  by  Equation  20.  If  the  width  of 
the  depletion  layer  is  related  to  the  voltage  at  the  source  edge 
of  the  depletion  layer  by 

(1  -  b2)  =  K/V(Xj)  ,  (22) 

then  the  current-voltage  characteristic  is  as  shown  in  Figure  6.  As 
the  channel  depletion  layer  gets  wider  with  increasing  drain  to  source 
voltage,  a  larger  fraction  of  the  total  voltage  drop  occurs  under  the 
depletion  layer  and  the  current  saturates.  The  channel  under  the  de¬ 
pletion  layer  does  not  pinch  off  in  this  approximation.  In  fact, 

B  pinchoff  or  bj  +  0  would  imply  an  infinite  potential  drop  between 
drain  and  source.  The  limitation  of  the  zero  space-charge  approxima¬ 
tion  can  be  seen  by  referring  to  Equation  20.  There  will  be  current 
continuity  in  the  channel  as  long  as  the  current  is  linearly  related 
to  the  field  or  the  mobility  is  constant.  In  Si  devices  the  mobility 
decreases  with  increasing  electric  field  as  follows: 

p(E)  =  pQ  exp  (-  y0E/vsat)  .  (23) 

At  low  fields  the  mobility  is  approximately  constant  and  there  is  no 
space  charge  in  the  channel.  As  the  drain-source  voltage  and  field 
in  the  channel  increase,  the  mobility  decreases  and  there  is  a  build 
up  of  electrons  to  satisfy  the  current  flow.  These  extra  electrons 
further  increase  the  potential  drop  in  the  channel  by  further  increas¬ 
ing  the  magnitude  of  the  electric  field  in  the  channel.  Beyond  the 
end  of  the  depletion  layer  toward  the  drain  end  of  the  channel  the 
channel  width  increases.  Since  the  field  and  total  current  must  be 
continuous  through  this  region,  the  electron  concentration  must  drop 
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FIG.  6  DRAIN-CURRENT  VS.  DRAIN-SOURCE  VOLTAGE  FOR  ZERO  SPACE-CHARGE  APPROXIMATION. 


below  the  background.  This  depletion  region  becomes  the  second  half 


of  the  channel  charge  dipole.  The  depletion  layer  lowers  the  magni¬ 
tude  of  the  channel  electric  field  which  increases  the  mobility  and 
reduces  the  amount  of  depletion  further  along  the  channel. 

The  operation  of  a  GaAs  transistor  is  slightly  different.  For 
GaAs,  the  mobility  is  approximately  constant  up  to  the  peak  velocity 
field.  Beyond  the  peak  velocity,  the  velocity  decreases  with  in¬ 
creasing  field.  For  low  fields  in  GaAs  FETs  the  amount  of  depletion 

mm 

or  accumulation  is  low,  typically  less  than  10  times  the  background 
charge.  At  a  field  slightly  above  the  peak  velocity  field  the  electron 
velocity  starts  to  decrease.  This,  in  turn,  results  in  an  accumulation 
formation  at  this  point  as  the  higher  velocity  electrons  from  the 
source  end  of  the  channel  arrive  at  the  high  field  point  and  slow  down. 
This  further  increases  the  field.  The  field  at  equilibrium  outside  the 
dipole  must  be  low  enough  to  prevent  further  buildup  of  charge.  This 
results  in  a  lowering  of  the  drain  current  after  the  dipole  has  formed. 
The  field  under  the  gate  depletion  layer  on  the  source  end  must  be 
lower  than  the  peak  velocity  field  to  prevent  a  second  dipole  from 
forming.  By  changing  the  drain-source  voltage,  the  magnitude  of  the 
field  under  the  entire  gate  can  be  pulled  above  tbe  peak  field.  How¬ 
ever,  the  electrons  will  still  accumulate  at  the  highest  field  or 
lowest  velocity  point  reforming  a  larger  dipole  at  this  point  and  again 
reducing  the  field  in  the  channel  below  the  peak  field.  Beyond  the 
narrowest  part  of  the  channel  the  magnitude  of  the  field  starts  to 
decrease.  This,  in  turn,  increases  the  electron  velocity  and  reduces 
the  electron  concentration,  forming  a  depletion  layer.  This  further 
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reduces  the  field  below  the  peak  field,  moving  the  electron 
concentration  back  to  the  doping  level  near  the  drain  contact. 

Comparison  with  Two-Dimensional  GaAs  Results.  The  recent  work 
of  Laux1 0  is  the  first  two-dimensional  simulation  of  a  GaAs  FET  with 
reasonable  doping  levels.  The  device  modeled  has  a  1-ym  long  channel, 
a  0.25-ym  long  gate,  and  a  0.25-ym  high  channel.  The  velocity-field 
curve  used  is  Equation  10  with  y  =  5000  cm2/V-s,  v  .  =  1  x  107  cm/s 

O  S&u 

and  E^  =  2690  V/cm.  The  diffusion-field  curve  used  is  Equation  11 
with  D  =20  cm2/s  and  E.  =  1.33  E  .  Figure  7  shows  the  electron 
concentration  and  potential  along  the  channel  for  the  two-dimensional 
model.  The  break  in  the  1.5-V  electron  concentration  is  numerical. 
Figures  8  through  10  show  similar  results  for  the  present  model. 

Figures  8  through  10  have  the  0. 2-,  0. 4-  and  1 . 5-V  information  for  the 
simple-parallel  and  sloping-potential  depletion-layer  approximations. 
The  geometry -distance  plots  in  the  three  figures  show  the  two  approxi¬ 
mations.  The  sloping-potential  depletion  layer  is  always  smaller  than 
the  simple-parallel  model.  It  is  also  sloping  upward  where  the  simple 
one  is  narrowest.  Figure  8  with  a  0.2-V  drain-source  bias  shows  the 
FET  at  its  peak  current  point,  with  the  velocity  at  the  peak  velocity 
of  1.2  x  107  cm/s.  The  velocity  at  the  drain  and  source  contacts  is 
approximately  the  same  for  the  two  depletion-layer  models.  For  all  the 
GaAs  devices  modeled,  the  choice  of  the  depletion-layer  approximation 
changed  the  distribution  of  electrons  and  field  in  the  channel  but  did 
not  greatly  change  the  terminal  drain  current.  Because  the  field  in 
the  channel  for  the  0.2-V  case  is  below  the  peak  velocity  field,  the 
electron  concentration  is  almost  equal  to  the  background  doping.  For 
all  the  cases  modeled,  the  simple  model  slightly  underestimates  the 
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FIG.  7  (a)  ELECTRON  CONCENTRATION  (m"’)  AND  (b)  POTENTIAL  (V)  FROM 

S’  TO  D'  ALONG  THE  y  =  0  SYMMETRY  LINE  FOR  THE  UNIFORMLY  DOPED 

INTRINSIC  MESFET.  V_,c,  =  0  V;  V„,_,  =  0.2,  O.k  AND  1.5  V, 
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field  in  the  channel  under  the  gate  and  overestimates  the  field  in  the 


rest  of  the  channel.  This  means  that  a  better  approximation  is  needed 
for  the  depletion-layer  geometry. 

Figure  9  shows  the  results  for  the  O.^-V  drain-source  bias 
condition.  Here  a  charge  dipole  has  formed  in  the  channel.  The  peak 

a 

electron  concentration  is  1.1  times  the  background  doping  and  the  mini¬ 
mum  concentration  is  0.7  times  the  background  concentration.  There  is 
good  agreement  between  the  simple-parallel  and  the  two-dimensional 
results.  The  peak  electric  field  is  approximately  20  kV/cm  and  occurs 
beyond  the  drain  end  of  the  metal  gate  as  shown  in  the  field-distance 
plot  of  the  figure.  The  excess  charge  region  extends  under  approximately 
two  thirds  of  the  gate. 

Figure  10  shows  the  results  for  the  1.5-V  comparison.  This 
condition  shows  the  effect  of  the  two  depletion-layer  approximations 
on  the  electron  concentration  along  the  channel.  The  simple  depletion- 
layer  approximation  gives  a  narrower  channel.  This  results  in  a  high 
electron  concentration  that  is  more  peaked.  The  sloping-potential 
approximation  gives  less  change  in  channel  geometry  with  position.  This 
spreads  out  the  excess  electron  peak  and  reduces  its  maximum  value. 

This  is  closer  to  the  two-dimensional  result  of  Figure  7.  The  sloping 
potential  also  lowers  the  peak  electric  field.  The  peak  field  has 
moved  closer  to  the  drain  contact  under  this  bias  condition.  The  field 
is  above  the  peak-velocity  field  through  the  entire  gate. 

Figure  11  shows  the  potential  vs.  distance  results  for  the  six 
cases.  The  sloping  potential  depletion  layer  gives  a  slower  rise  in 
potential  and  moves  the  curve  to  the  left.  Again,  these  curves  agree 
well  with  the  two-dimensional  results  of  Figure  7b. 
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The  common  source  output  characteristics  of  the  two  models  are 
shown  in  Figures  12  and  13.  The  simple  model  drain  currents  are 
approximately  25  percent  higher  than  the  two-dimensional  results.  This 
difference  is  a  result  of  the  high  mobility  of  GaAs.  The  600-mA 
difference  between  the  drain  currents  at  zero  gate  bias  corresponds  to  an 
electric  field  difference  of  only  300  V/cm  at  the  contacts.  This  is 
less  than  1  percent  of  the  peak  channel  field.  Further  study  of  the 
depletion-layer  model  and  the  two-dimensional  results  are  needed  to 
resolve  this  difference.  Although  the  drain  currents  in  the  simple 
model  are  high,  the  change  with  gate  bias  is  approximately  correct. 

The  drain  current  for  a  0.5-V  drain-source  voltage  at  different  gate 
voltages  is  shown  in  Figure  14.  Over  a  gate  voltage  range  from  zero 
to  -2  V  the  slopes  of  the  curves  and  thus  the  transconductance  are 
approximately  equal.  A  similar  comparison  is  shown  in  Figure  15  for 
a  1.0-V  drain-source  voltage. 

Summary  of  Comparison.  Over  the  range  of  data  available  there 
is  good  agreement  between  the  present  model  and  a  full  two-dimensional 
simulation.  There  is  very  good  agreement  on  the  potential  and  electron 
distribution.  Although  the  currents  in  the  simple  model  are  high,  the 
electric  field  in  the  channel  is  only  off  by  300  V/cm.  The  best 
approach  at  this  time  would  be  to  make  further  runs  with  both  programs 
in  order  to  pinpoint  the  reasons  for  the  differences.  However,  there 
are  two  problems  with  this.  First,  the  two-dimensional  problem  has 
a  relatively  large  mesh  spacing.  It  is  not  clear  that  further  runs 
would  give  much  additional  information.  Second,  the  two-dimensional 
program  is  costly.  A  current  voltage  set  such  as  shown  in  Figure  12 
costs  hundreds  of  dollars  of  computer  time.  By  comparison,  the  data 
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FIG.  15  TRANSCONDUCTANCE  COMPARISON. 


in  Figure  13  cost  less  than  $3.  It  would  be  better  to  use  the  simple 
model  data  as  input  to  the  full  program.  Finally,  although  the 
currents  are  high,  the  transconductance  is  a  more  important  parameter 
for  ac  operation. 
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SECTION  k 


EFFECTS  OF  VARIOUS  PARAMETERS  ON  FET  OPERATION 

Introduction.  The  idea  behind  the  study  of  FETs  is  to  understand 
better  their  operation  so  that  their  actual  performance  can  be  improved. 
Subject  to  technical  limitations,  the  device  builder  has  a  choice  of 
materials  and  geometry  when  building  FETs.  The  material  choice,  in  turn, 
gives  a  particular  velocity-field  and  diffusion-field  characteristic,  in 
this  section  a  variety  of  geometries,  velocity-field  choices,  and  diffusion- 
field  choices  are  investigated  in  order  to  provide  a  better  understanding 
of  FET  operation. 

Gate  Length  Variation.  The  effect  of  gate  length  on  FET  operation 
is  considered  here.  The  device  is  1  ym  long  with  a  0.25-ym  high  channel 
doped  at  lxl0I7/cm3.  The  velocity-field  and  diffusion-field  charac¬ 
teristics  are  given  by  Equations  10  and  11  with  yQ  =  5000  cm*/V-s, 

v  .  =  107  cm/s,  E  =  3500  V/cra,  E.  =  1.33  E  .  The  resulting  charac- 
sat  ’v  d  v 

teristics  for  gate  lengths  of  0.01,  0.2,  0.3  and  0.5  ym  are  shown  in  Figures 
l6  through  19.  The  drain  current  vs.  drain-source  voltage  at  zero 

gate  bias  is  shown  in  Figure  20  for  the  various  cases.  At  channel 
fields  below  the  peak  velocity  field  there  is  little  space  charge  and 
Equations  l8  and  19  apply.  The  current  depends  on  the  boundary 
electric  field,  so  increasing  the  gate  length  lowers  the  current. 

However,  even  at  a  very  short  gate  length  of  0.01  ym,  the  gate  deple¬ 
tion  layer  covers  20  percent  of  the  channel  length.  Since  the  integral 
of  the  total  field  along  the  channel  in  each  case  must  equal  the 
drain-source  voltage,  the  narrower  gate  gives  a  higher  field  and 
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FIG.  1 6  CHARACTERISTICS  OF  O.Ol-Mm  GATE  FET. 
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FIG.  IT  CHARACTERISTICS  OF  0.2-pm  GATE  FET. 


DRAIN- SOURCE  VOLTAGE 


FIG.  20  DRAIN  CURRENT  VS.  DRAIN-SOURCE  VOLTAGE  WITH 
GATE  LENGTH  AS  PARAMETER. 
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current  at  the  peak  current  point  around  a  drain-source  voltage  of  0.2  V. 
Increasing  the  gate  length  decreases  the  peak  current.  Beyond  the 
peak  current  point  a  charge  dipole  forms  for  each  gate  length.  After  the 
charge  dipole  forms  it  dominates  the  current  flow  in  the  channel.  Ex¬ 
cept  for  the  position  along  the  channel,  the  dipole  in  each  case  is 
approximately  the  same.  The  peak  field  in  the  dipoles  for0.5-V  drain- 
source  voltage  varies  less  than  5  percent  for  the  range  of  gate  lengths 
considered.  The  maximum  and  minimum  charge  concentrations  vary  even  less. 
Beyond  0.3-V  drain-source  voltage  on  Figure  20,  the  drain  current  is 
approximately  the  same  for  each  case.  Since  the  current  characteris¬ 
tics  of  FETs  with  varying  gate  lengths  are  the  same,  the  main  effect  of 
varying  the  gate  length  is  the  capacitance  of  the  gate.  Changing  the 
gate  length  changes  the  volume  of  the  depletion  layer  and  this  changes 
the  amount  of  current  required  to  change  the  layer.  However  even  the 
very  short  gate  0.01-pm  device  occupies  approximately  20  percent  of 
the  channel  length  at  zero  gate  bias. 

Diffusion-Coefficient  Variation.  Figures  21  through  25  show  a 
0.25-pm  gate  FET  with  the  same  velocity-field  characteristic  and  geometry 
as  the  gate  length  variation  test.  The  diffusion  coefficient  vs.  field 
curve  has  been  multiplied  by  a  factor  of  0.5,  1,  5,  10  and  20  in  the 
figures.  Figure  2 6  shows  the  resulting  drain  current  vs.  drain-source 
voltage  curve  for  the  various  diffusion  coefficients.  Progressing 
through  Figures  21  through  25  shows  that  the  effect  of  increasing  the 
diffusion  coefficient  is  to  reduce  the  gradient  in  the  electron  con¬ 
centration.  This  lowers  the  magnitude  of  the  maximum  and  minimum 
electron  concentrations,  spreads  out  the  charge  dipole,  and  spreads 
out  and  lowers  the  electric  field.  For  drain-source  voltages  lower 
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FIG.  22  DIFFUSION  VARIATION  TEST  FACTOR 
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FIG.  25  DIFFUSION  VARIATION  TEST  FACTOR  =  20. 
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than  0.25  V,  the  electron  concentration  is  approximately  equal  to  the 
background  doping,  so  there  is  no  electron  gradient  and  thus  no  diffusion 
current.  The  low  voltage  current  is  therefore  the  same  regardless  of 
the  diffusion  coefficient.  Above  0.25-V  drain-source  voltage,  the 
current  depends  slightly  on  the  diffusion  coefficient.  The  electric 
field  curves  shown  in  Figure  21  for  the  0.5  factor  case  have  a  lower 
potential  drop  across  the  charge  dipole  than  the  rest  of  the  curves. 

This,  in  turn,  increases  the  field  outside  the  dipole  and  increases 
the  current.  The  factor  =  1  and  factor  =  5  curves  are  also  plotted  in 
Figure  26.  They  are  approximately  the  same.  The  factor  *  10  and 
factor  =  20  curves  are  the  same  as  the  factor  =  5  curve. 

From  this  set  of  data  it  is  found  that  diffusion  coefficient  mag¬ 
nitudes  greater  than  the  typical  value  by  factors  up  to  20  have  no 
effect  on  the  terminal  current  voltage  characteristic.  The  potential 
drop  over  the  resulting  dipole  is  the  same  even  though  the  shape  and 
magnitude  of  the  field  and  electron  concentrations  change.  Thus,  effects 
which  increase  the  diffusion  coefficient  but  do  not  change  the  velocity- 
field  curve  should  not  affect  the  FET  terminal  I-V  characteristic. 

Velocity-Field  Curve  Variation.-  In  this  section  the  effect  of 
the  velocity-field  characteristic  on  FET  operation  is  studied.  The 
velocity-field  curve  used  is  Equation  10: 


v(E) 


1*„E  *  vsat  (E/E»)“ 
1  +  IE/Ev)l» 


For  vsa1.  =  107  cm/s  and  Ey  =  3000  V/cm,  the  velocity-field  char¬ 
acteristics  for  mobilities  of  1*000,  5000  and  6000  cm2/V-s  are  shown  in 
Figure  27.  The  resulting  drain  current  vs.  drain-source  voltage  char¬ 
acteristic  is  shown  in  Figure  28.  For  bias  voltages  below  0.2  V, 


VELOCITY,  10  cm/s 


FIG.  27  VELOCITY-FIELD  CHARACTERISTICS  FOR 
VARIOUS  MOBILITIES. 
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FIG.  28  FET  OPERATION  WITH  VARYING  MOBILITY. 
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there  is  no  space  charge  in  the  channel  so  the  electric  field  in  the 


channel  depends  only  on  the  geometry  and  is  the  same  for  the  three 
mobility  cases.  The  current  difference  then  depends  only  on  the  mobility. 
Above  the  bias  for  peak  current  (0.2  to  0.25)  the  higher  mobility  currents 
are  greater  than  the  low  mobility  ones  but  the  variation  is  much  less 
than  the  ±20-percent  variation  in  the  low  field  mobility. 

Saturat ed  Velocity  Variat ion.  For  =  5000  cm2/V-s  and  E^  =  3000 
V/cm,  the  velocity-field  curves  for  v  =  8  x  106,  107  and  1.2  x  107 

S3,  u 

cm/s  are  shown  in  Figure  29.  The  corresponding  current-voltage  charac¬ 
teristic  is  shown  in  Figure  30.  Here,  since  the  low-field  velocity-field 
curve  is  the  same  in  each  case,  the  linear  portion  of  the  curve  is  also 
the  same.  Again,  a  + 20-percent  change  in  the  saturated  velocity  produces 
less  than  a  10-percent  change  in  the  drain  current  in  the  saturated  region 
of  operation. 

Peak  Field  Variation.  For  u  =  5000  cm2/V-s  and  v  ^  =  107  cm/s, 

-  "  O  SSL  U 

the  E^  was  varied  between  2000  and  5000  V/cm.  The  velocity-field  curves 
for  these  cases  are  shown  in  Figure  31  and  the  resulting  drain  current 
vs.  drain-source  voltage  curves  are  shown  in  Figure  32.  Here,  both  the 
low-field  and  high-field  velocities  are  the  same.  However,  there  is  a 
greater  than  ±10-percent  variation  in  the  resulting  current.  This  leads 
to  the  conclusion  that  the  saturation  current  depends  not  on  the  low  field 
mobility  or  saturated  velocity  directly  but  on  the  peak  velocity.  To 
show  this,  the  current  for  a  drain-source  voltage  =  0.35  V  was  plotted 
vs.  the  peak  velocity  for  the  various  mobilities,  saturated  velocities, 
and  peak  fields  in  Figures  28,  30  and  32.  This  is  shown  in  Figure  33. 

This  i'.iows  the  linear  relationship  between  peak  velocity  and  saturated 
current.  That  the  saturated  current  should  depend  on  the  peak  electron 
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FIG.  29  VELOCITY -FIELD  CHARACTERISTICS  FOR 
VARIOUS  SATURATED  VELOCITIES. 
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FIG.  30  SATURATED  VELOCITY  VARIATION  EFFECTS  IN  FETs 


FIG.  31  VELOCITY-FIELD  CHARACTERISTICS  WITH  E  AS  A  PARAMETER. 
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FIG.  32  PEAK  VELOCITY  VARIATION  EFFECTS  IN  FETs. 
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FIG.  33  SATURATED  CURRENT  VS.  PEAK  VELOCITY  IN  FETs 


velocity  rather  than  on  the  mobility  or  the  saturated  velocity  is 
reasonable  for  GaAs  FETs.  For  fields  below  the  peak  velocity,  there  is 
little  space  charge  in  the  channel.  The  high  peak  velocity  in  the  high- 
field  part  of  the  channel  corresponds  to  the  high  velocities  and  thus 
high  currents  at  the  drain  and  source  contacts.  Beyond  the  peak  velocity 
field  point,  the  current  begins  to  drop  with  increasing  bias.  The 
amount  of  current  dropback  depends  on  the  velocity  dropback.  This  is 
shown  in  Figure  3^  for  the  various  cases  of  mobility,  saturated  velocity, 
and  peak  field.  The  larger  velocity  dropbacks  give  larger  peak  to 
saturation  current  dropbacks,  but  the  larger  peak  velocities  result  in 
larger  saturated  currents. 

Summary  of  Parameter  Var iati'-'n.  FET  operation  over  a  range  of 
gate  lengths,  diffusion  coefficients,  and  velocity-field  characteristics 
have  been  studied.  The  main  factor  in  the  saturated  current  and  thus 
the  transconductance  is  the  peak  velocity  of  the  material  used.  The 
gate  length  and  diffusion  coefficient  have  almost  no  effect.  The  mobili¬ 
ty,  saturated  velocity,  and  peak  velocity  affect  the  current  by  affect¬ 
ing  the  value  of  the  peak  velocity.  Although  the  current  is  not  affected 
by  the  low  field  mobility,  the  series  resistance  and  thus  the  parasitics 
will  be.  Also  the  gate  length  affects  the  gate  capacitance  and  thus  the 
RF  operation  of  the  FET. 
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SUMMARY  AND  CONCLUSIONS 

The  modeling  and  operation  of  GaAs  FETs  has  been  discussed  in  this 
report.  After  discussing  the  basis  of  the  model,  the  results  were  compared 
with  similar  two-dimensional  results.  There  is  very  good  agreement  between 
the  electric  field  and  electron  distributions  in  the  two  models.  The 
drain  currents  for  the  simple  model  are  approximately  25  percent  higher 
than  the  similar  two-dimensional  results.  This  problem  can  be  solved 
with  an  improved  depletion-layer  model. 

Various  operating  characteristics  of  FETs  were  discussed  in  Sec¬ 
tion  k.  The  effect  of  velocity-field  and  diffusion-field  characteris¬ 
tics  on  device  operation  were  studied.  The  diffusion  coefficient  has 
little  effect  on  the  current-voltage  characteristic  of  the  device.  The 
main  parameter  affecting  the  current  in  the  saturated  region  of  operation 
is  the  peak  electron  velocity.  The  gate  length  has  little  effect  on  the 
current  in  the  saturated  region.  This  is  because  the  current  flow  is 
dominated  by  the  charge  dipole  formed  for  saturated  operation.  The 
simple  model  developed  here  can  be  extended  to  study  the  following 
important  parameters  concerning  FET  operation: 

1.  Small-signal  characteristics. 

2.  Electric  field  breakdown. 

3.  High-power  design. 

i*.  Nonequilibrium  electron  transport  in  the  Channel. 

5.  Large-signal  operation. 

6.  Noise  performance. 
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This  can  be  achieved  as  follows.  After  the  dc  solution  for  a  given 
FET  geometry,  material,  and  operating  point  has  been  found,  the  small- 
signal  operating  conditions  as  a  function  of  frequency  can  also  be 
determined.  This  involves  making  a  linear  expansion  about  the  dc  solution 
and  solving  the  resulting  matrix  equation.  The  small-signal  equations  have 
been  written  as  a  subroutine  for  the  dc  program.  The  program  is  pres¬ 
ently  being  added  to  the  dc  program  and  should  be  available  soon. 

With  the  dc  solutions  available,  the  conditions  of  operation  leading  to 
avalanche  breakdown  can  be  studied.  There  are  two  high-field  regions 
in  the  device,  at  the  drain  end  of  the  gate  metal  and  in  the  narrowest 
part  of  the  conducting  channel.  The  ionization  integral  can  be  found 
under  various  bias  conditions  to  find  the  limits  of  high-power  operation. 
This  is  an  extension  of  the  electric  field  breakdown  problem.  The  current 
in  the  channel  depends  on  the  material  parameters  and  geometry  as  discussed 
in  this  report.  A  variety  of  other  possible  FET  materials  and  struc¬ 
tures  needs  to  be  studied  to  find  the  best  high-power  structure. 

The  area  of  hot-electron  effects  can  also  be  investigated.  For  a 
rapidly  changing  field  (both  in  space  and  time)  in  the  channel,  the 
electron  velocity  is  not  well  known.  This,  in  turn,  will  affect  the 
current  in  the  channel  and  the  operation  of  the  FET. 

The  FET  program  as  written  is  a  time-dependent  solution.  The 
dc  solution  is  found  by  solving  the  current  transport  and  Poisson's 
equations  through  a  transient  to  a  steady-state  solution.  As  presently 
written,  the  depletion-layer  charging  current  is  not  included.  By 
adding  this  term  to  the  equations,  large-signal  operation  and  modeling 
can  be  achieved. 
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Again,  as  in  the  case  of  small-signal  operation,  having  a  dc 
solution  makes  a  small-signal  noise  model  possible.  This  model  could 
improve  on  the  simple  analytic  noise  models  presently  available  because 
of  the  increased  complexity  of  the  underlying  dc  solution. 
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